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Abstract 

o 

We give a definition and derive the equations of motion for the center of mass and angular momentum 
of an axially symmetric, isolated system that emits gravitational and electromagnetic radiation. A central 
{^j^ feature of this formulation is the use of Newman-Unti cuts at null infinity that are generated by worldlines 

of the spacetime. We analyze some consequences of the results and comment on the generalization of 
this work to general asymptotically flat spacetimes. 

1 Introduction 

The notion of center of mass for an isolated system is very important in newtonian theory. It is used to 
define the finear and intrinsic angular momentum of the system, both conserved observables in the theory. 
However, its generalization to General Relativity (GR) has proved to be a non trivial task. 

A major obstacle to a relativistic definition is the fact that energy or momentum cannot be local quan- 
tities. One cannot use the Energy Momentum tensor of a system as in other theories since gravitational 
waves carry away energy and momentum and nevertheless are solutions of Ricci flat equations. Thus, one 
must search for global definitions of these quantities. Also one must bear in mind that, unless the spacetime 
is stationary, energy or momentum of an isolated system are not conserved in GR due to the emission of 
gravitational radiation. 

Fortunately there are ways to overcome, at least in principle, these difficulties. One has available in 
the literature the notion of an asymptotically flat spacetime, the appropriate framework to analyze isolated 
systems in GR. For those spacetimes one can define the notion of Bondi mass and linear momentum and 
write down equations of motion linking their time evolution with the emitted gravitational radiation. All 
that is needed then is to relate these variables to a suitable definition of center of mass, characterized by a 
worldline such that = MR°'+ radiative corrections. The problem is how to select this worldline. 
There is available in the literature two nice approaches based on information that can be retrieved at null 
infinity J^. 

One approach that has been carried out by Newman and collaborators is based on the introduction of 
asymptotically shear free null congruences, i.e., congruences such that at have vanishing shear. At null 
infinity this congruence appears to come from a point in Minkowski space. This asymptotically shear free 
condition yields a family of two-surfaces, "good cuts", that are constructed from special solutions of the 
Good Cut equation. This family is characterized by complex worldlines in a fiducial holographic space. 
Furthermore, the vanishing of the complex mass dipole term at singles out a particular worldline of this 
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family. Its real part gives the center of mass while the imaginary part is by definition the intrinsic angular 
momentum per unit mass. Assuming a quadrupole radiation and using any available definition of total 
angular momentum (where all of them coincide) one obtains equations of motion coupling center of mass, 
intrinsic angular momentum and radiation. (A complete description is available at a Living Reviews [T]). 

The other approach has been done by Moreschi [[2|. In this case one first defines the notion of "super- 
momentum" at J" and then asks for a special family of cuts, called nice cuts, where the supermomentum 
only has an / = and I = 1 spherical harmonics decomposition. Again one obtains another holographic 
solution space and one special worldline in this space is selected, via a similar condition as above, as the 
center of mass. Furthermore. Moreschi defines the notion of total angular momentum |l3| and its restriction 
to the center of mass worldline yields the intrinsic angular momentum of the system. 

Both formulations coincide at a linear level if the gravitational radiation is pure quadrupole. In spite 
of the clever ideas used in both approaches to define a global notion of center of mass there are some 
drawbacks that one must mention since in the end the definition should be used in gravitational wave 
astronomy to characterize the dynamics of compact highly energetic objects like AGN or binary systems 
such as BH-BH or BH-NS. 

It is not true that light coming from distant isolated systems is asymptotically shear free. In fact the 
shear of light coming from these sources is used to define weak lensing effects in GR. Most important, 
neither the nice cuts nor the good cuts are related to the future light cones from points inside the spacetime. 
One can write down the equation that must be satisfied by any light cone cut at coming from a worldline 
in the spacetime Q. None of these approaches satisfy this equation up to second and higher orders. 

The idea of this note is to use the light cone cut equation recently obtained m together with ideas 
borrowed from the two approaches to define center of mass and intrinsic angular momentum. So far, there 
is not available in the literature a satisfactory definition of angular momentum for non-stationary spacetimes 
without symmetries. To overcome this difficulty we will consider here axially symmetric spacetimes since 
in this case there is a suitable definition of angular momentum. 

This work is divided in five sections and two appendices. Section 2 is devoted to the mathematical tools 
and definitions needed for this work. Readers familiar with the Newman-Penrose formulation, asymptotic 
flatness and the Komar integral may skip this section. Sections 3 and 4 constitute the core of the paper 
Finally, in the Conclusions we summarize our results and outline a generalization of this formalism to 
arbitrary spacetimes. 

2 Foundations 

There are many result that are needed for this work. In this section, we introduce several of the key ideas 
and the basic tools that are indispensable in our later discussions. A thorough derivations of these results 
are given in the references. 

2.1 Asymptotically flat spacetime and 

The notion of asymptotically flatness is the adequate tool to analyze the gravitational and electromagnetic 
radiation coming from an arbitrary compact source. A spacetime can be thought of as asymptotically flat 
if the curvature tensor vanishes as infinity is approached along the future-directed null geodesies of the 
spacetime. All the null geodesies and up at what is referred to as future null infinity, the future 

boundary of the spacetime |[T], ||5]- We introduce a natural set of coordinates in the neighborhood of 
called Bondi coordinates {ub, r, (). In this system, the Bondi time ub labels a special family of null 
surfaces whose intersection with are two spheres, r is the affine parameter along each null geodesic of 
the constant ub surface and ( ~ e"^ cot |, is the complex stereographic angle that labels the null geodesies 
of the null surface. 

Associated with the Bondi coordinates is a null tetrad system based on these outgoing null hypersurfaces 
labeled by (la, Ua, rua, rria)- The first tetrad vector la is defined as Q 

L-VaUB. (2.1) 

Thus, 1°- = g°-^VbUB is a null vector tangent to the geodesies of the surface. For the second tetrad vector 
we pick a null vector n° normalized to 1°^ 

rial" = 1 (2.2) 



The tetrad is finally completed with the choice of a complex null vector to" orthogonal to /" and n° 

mam'' = -1 (2.3) 

The spacetime metric is then Q 

gab = Lnt + nJb - mamb - mamb (2.4) 

There is a great number of tetrad freedom, but the most important for us is a different choice of the original 
ub ~ const, cuts of ^+ so that 

UB = Z{u,CX) (2.5) 
where Z{u, (, () is a real function. Let us denote by T the inverse function Z, so 

u^T{uB,C,C) (2.6) 

is easy to show that T = then the rest of the coordinate system and the tetrad system are then constructed 
as before. 

A second freedom introduces the concept of spin weight Q. A quantity rj that transforms as 77 — > e^^^rj 
under a rotation to" c^^to" is said to have a spin weight s. For any function /(w, C,, Q, we can define 
two differential operators 3 and 3 by 

3/ = P'-'%^ (2.7) 

g/ = p^+-^^i£ZlIl (2.8) 

dC 

where / has a spin weight s and P is the conformal factor defining the sphere metric 

note that for axial symmetry, the operators 3 and 3 act as derived on 9. In addition, all functions will be 
functions that do not depend of (f>. In Eqs. (12.7b and (12.8b . the conformal factor P is arbitrary. However, in 
Bondi coordinates, the conformal factor is restricted to 

P = Po = l + CC (2.10) 

2.2 The Newman-Penrose formalism 

Although the Newman-Penrose (NP) formalism is the basic working tool for our analysis, we will simply 
give an outline of the formulation and leave the reference Q for details. We focus in the general form of 
the asymptotically flat solutions of Einstein-Maxwell equations in Bondi coordinates. 

The NP version ||5l,|l6l of the vacuum Einstein (or the Einstein-Maxwell) equations uses the tetrad 
components 

A"c = (^°,n°,m",TO"); c= 1,2,3,4 (2.11) 
rather than the metric, as the basic variable. The metric, Eq. ( 12.4b can be written as 

g-b^jfdx-^X\ (2.12) 

with 
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The Ricci rotation coefficients 7"^^^ are defined by ||5l,|l6l 



(2.13) 



(2.14) 



so 

Icdf = -Idcf (2.15) 
the 12 spin coefficients are defined as combinations of the 7*^^/: 
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The third basic variable in the NP formalism is the Weyl tensor or, equivalently, the following five 
complex tetrad components of the Weyl tensor: 

02 = -\{Cai,al''nH'n'' - Cabcdl'^n^i'm'') (2.17) 
03 = Cabcdl^n^n^m:^; 04 = -Cabcdn''m!'n''rn'^ 

When an electromagnetic field is present, we include the complex tetrad components of the Maxwell 
field 

00 = Fabl'^m'- 01 = ^Fabil'^n'' + 02 = Fabn''m\ (2.18) 

into the equations |[l], IS). 

The Peeling theorem of Sachs Q tell us that the behavior of the Weyl scalar and the maxwell tensor is 
given by: 



-00 




-5 _ 






^1 


= 0?r" 


-4 _ 


hO(r- 




-02 


= 02"r- 


-3 _ 


hO(r- 


-4) 


-03 


= 03^" 


-2 _ 


hO(r- 




V'4 


= ^Ir- 


-1 _ 


hO(r- 




00 


= 0°r- 


-3 


hO(r- 




01 


= 0?r- 


-2 _ 


hO(r- 




02 


= 0^^- 


-1 







where the quantities with a zero superscript are function only of {ub, C, 0- The spin coefficients and metric 
variables are given as IH, Q. 

K = n = e = 0; P^p; T~a + l3 

p = -r-^ - a^a^r-^ + Oir--') 

a = (7"r-' + [((70)2^0 - ^;j/2]r-4 + 0(r-5) 

a = a°r-i +C'(r-2) 

/3 = /3°r-i +0(r-2) (2.20) 

7 = 7" -?A^(2r2)-i +0(r-3) 

A = A"r-i+0(r-2) 



where the relationships among the r-independent functions 
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(2.21) 

Finally we return to the choice of the null tetrad. If we start from the rescaled metric g"''* = Z'^ g"'' the 
tangent vector to the generators of is rescaled as — Z'n°^ . Using this null vector we can define all 
the other vectors of the new tetrad as 

= -^{la- —ma- —ma + ^Ua) (2.22) 

Z' r r 

K = -^na (2.23) 

w* = -TTii'ma- —ria) (2.24) 
Z r 

= -^.{ma- -ria) (2.25) 
Z' r 

t3 T 

L{ubX,0 = " ^iu)Z{u,C,C)\u=TiuBX,0- 

In the above 3( [,a) or g(-„) means to apply the eth operator keeping ub or u constant respectively. 
From this tetrad we can define the new Weyl scalar Ul, ifTTI 

<KC,C) ^ [^,o_3^^o^3^2^o_^3^0](^^^^^^) (2.26) 



where 



Z'2 

from which we can write the following approximation 



V'?* - ^? - 3X^2 (2.27) 
if we keep up to linear terms in Z' and/or L. 

2.3 Evolution equations and physical definitions 

Using the peeling theorem all the radial part of Einstein equations can be integrated leaving only the 
Bianchi identities at as the unsolved equations. Some of those equations are used to relate Weyl scalars 
with the free Bondi data, i.e., fS] 

V'^-V^ = gV"-g2^" + a"a"-a"^° (2.28) 
V'3 = 3^° (2.29) 



^2 = (2.30) 
In the above a'^ is the free data. From Eq. (12.28b we can define the so called mass aspect HI 

* = V2° + + f^"^ (2.31) 

which satisfies the following condition 

^'=W (2.32) 
and finally the evolution equations (Bianchi identities) Q 

vi? = -g* + g-V + ga"^° + 3a"g5^° + ^0?02 (2.33) 



2G ,n-0 



= -g^g^" - cr°a ' + ^0^02 (2.34) 



1^2 — —uu—uu^ — V'2V^2 

6? = -g<^^ 



(2.35) 

00 = -d4 + a"4 (2.36) 



Using the mass aspect with -02' the second of the asymptotic Bianchi identities can be rewritten in the 
concise form 

* = -"^" + ^^2^- (2.37) 

Note that in the above equation the gravitational radiation cr" and the electromagnetic radiation determine 
the mass aspect ^P. In addition, we can define the Bondi mass and linear momentum as 

M . -5^/«0 (238, 

P' = -l^j m'dil, (239) 

and one can easily see that the Bondi mass decreases as a result of the emitted radiation. 



3 Angular Momentum 

The definition of angular momentum in general relativity has proven to be a major task which so far does 
not have a satisfactory solution. Basically the problem lies at identifying a canonical origin at null infinity. 
However, for vacuum axially symmetric spacetimes one can use the Komar integral associated with the 
rotation KilHng field and write a conserved quantity 

1 



r = — Mm <pV''eMdSab = const. (3.1) 

St 

We now want to extend this definition to include the contribution of the electromagnetic radiation. 
Using Stokes theorem and the fact that ^"^j is a Killing field we have 

j^^S,\^)dSab = 2 j RabS.\^)dY.'' (3.2) 
as ^ 

where 91] is the boundary of the hypersurface S. Since the Killing vector is tangent to ^ , ie ^^^^^ni, = 0, 
we can replace the Ricci tensor by the stress energy tensor Tab, in the above, i.e., 

j^'^S.UdSab = 167r j Tab^Md^" (3.3) 
as ^ 

Inserting the stress-energy tensor of electromagnetic field Tab ~ Tni^acF^ — \gabF^'^Fcd) in the r.h.s. of 
the above equation yields 

j^W^^i^^dSab = 4 J Fa''Fbc£.\^)dY.^' (3.4) 
as ^ 

= 4 / F/(VbA, - V ,Ab)£,t.dY.'' (3.5) 
Js 

Since we can choose the Maxwell potential to have axial symmetry and the Maxwell field is pure radiation 
we have in addition 

and 

V'Fac = 0. 

Thus, 

jv^'il^^dSab = -4y V Mbi\^^Fa')d^r (3.6) 
as ^ 



Using Stokes theorem once again, we finally obtain 



^4 J VMbiUF^d^a - -2(^A,^l^F''-dSac. (3.7) 
We thus redefine the angular momentum Eq. ( 13.1b . to include electromagnetic field, as 

^ islToo / [^"^(-) + ^MUF'^']dSa, (3.8) 

St 

The first integrand is the original gravitational term whereas the second one is the electromagnetic contri- 
bution. Using the N-P formalism||5l one can write Eq. ( 13.81 ) as (see Appendix A), 

J| = ^[I^i^i + - 4/m(AV?)|/=i] (3.9) 

where A'^{ub, Ci is the Maxwell potential free data related to the electromagnetic radiation via 

The new conserved quantity J|. will be called total angular momentum for any axially symmetric 
Einstein-Maxwell spacetime ie 

j| = (3.10) 

4 Center of mass 

Since by assumption the spacetime is axially symmetric we will also assume the center of mass is given 
by a worldline R"-{u) along the axis of symmetry, i.e., along the z-axis. We recall that in section l272l we 
introduced a null tetrad based on outgoing null hypersurfaces u = const. We will then assume that this 
family of hypersurfaces has been generated by the future light cones of R°-{u). The intersection of these 
light cones with ^+ yield Newman-Unti coordinates (w, C, Q. The basic idea is to start with the mass 
dipole term at in a Bondi frame, use Eq. ( 12.27b to write down the transformation equation to a Newman- 
Unti frame and demand that the mass dipole term vanishes on the u ~ const, slices. (A similar idea has 
been used before in the Kozameh-Newman approach for axially symmetric spacetimesfS] but as we will 
see later following both approaches yield different results.) 

4.1 Analysis and definition 

In a Bondi frame, the mass dipole momentum for asymptotically flat spacetime is defined to be the real 
part of the I = 1 component of i?e[i/;"]. We extend this definition to a Newman-Unti frame and define the 
mass dipole momentum as the I = 1 component of i?e[?/)]'*]. 

The basic idea to obtain the center of mass is to start by imposing the condition that on the u = const 
cuts generated by the worldline R"-{u), the mass dipole momentum vanishes. Then, using the relation (12.27b 
and expanding ij^l in a tensorial spherical harmonic basis as 

i;l = i,l\uB)Yl^ + rPl''{uB)Yl^+... 

one should obtain a relationship between ^[psi\^{uB)] and the center of mass worldline B}. 
Following this prescription, and using Eq. (l2.27b . on a u = const, slice we impose 

Re[i:l{Z, C, ' 3g(Z)^0(Z, C, OlUconst. = 0, (4.1) 

where we have replaced ub by the function ub = Z{R°'{u)XtC}- Furthermore, using a slow motion 
approximation and keeping up to first order terms in the velocity of the center of mass we write 



Z(i?"(u),C,C) ^u + 5u, 



where we assume 5u is small. Thus, we make a Taylor expansion of 

i?e[V'?*(M, C, C)] = Re[^l{u + 5u, (, C) - S5Sui^^^{u + Su, C, C)], 

decompose each term in spherical harmonics and demand that on the u ~ const, cut the I = 1 part of this 
series vanishes. The Taylor expansion yields 

Re[^'l*{u, C, C)] = Rei^Pliu, (, + V'?'(«, C, C)Su - 3c5Sui;°{u, C, C)] 

where we have omitted second order terms in Su. Taking the I = 1 part of the above expression and putting 
it equal to zero yields the following expression, 

Re[ij°^{u)Y = i?e[(g* - d^a°)5uY + 3i?e[35M(«' - d'^a°)]\ (4.2) 

i.e., the real, I = 1 part of V-'i can be written in terms of the center of mass and other Weyl scalars at null 
infinity. 

Inserting the following tensorial spin-s harmonics expansion lfT2ll 

Z = u + Su 

Su = -\r'{u)Y,°.{C) + ^''{u)Y,%{C) + x^^'{u)Y,%,{0 
dSu = {u)Y,\ (C) - 6x'^ (C) - I2x^'' iu)Yi,^, {() 

as = a'^ {ub)YI^ (C) + o^''' (ub)YI^^ {() 
^0 = i,f{uB)YlAO + i^i'{uB)Y^.,{C) + ^':'\ub)Y^,,AO 

* = -'^-^M{uB)-^-^P\uB)Y^AC)+-^''{uB)Y^,,{0 + -^''HuB)Yg,u{^^ 



V2 

in Eg. (14.21 1 yields 



I ^ 4>1\ub)yUC) + 4>1''{ub)yI,{0 

^ Q{ub) + ^f{uB)Y^m + 4>'u,{nB)YZ,{C) 
P - <Pf{uB)Y,-\0+c^f'iuB)Yi,l{0 



Qfi "^r^ 

D' = MR' ^x'^P^ + ^—Ri-qi'^ 

5V2c 5V2G 



7V2G 



{x'^{J2a%^ - 216^'^'^) + ^'^'=(432*^^ - 360(t^^')], (4.3) 



where we have defined 



(Since we assume axial symmetry, R' (u) only has a z component. Likewise, all the high order tensors are 
symmetric, diagonal and trace-free.) 

We now take a small digression to concentrate on the light cone cut function Z = Su defined as the 
intersection of the future lightcone from a worldline with The function Z dynamically depends 

on the matter and radiation content of the spacetime via the solution of the Einstein equations, i.e., the light 
cone cut function is dynamical variable and we do not make any a priori assumption about its behaviour. Z 
satisfies the equation 

s^z = A(z, gz, gz, ggz, c, C), 

and A satisfies the Einstein's equations (it vanishes for a flat spacetime). The freedom in the solution is 
given by a combination of Z 0, 1 spherical harmonics since they are annihilated by the g^ operator 
One can write this freedom as 

= t{u) + (C, C) - « + ^\u)Y^^. 



where int he last equality we have thrown away quadratic and higher order terms in The function 

A only contains I = 2 and higher spherical harmonics decomposition which are completely determined 
from the Einstein's equations. For example, if we assume a vacuum spacetime in the neighborhood of null 
infinity, the linearized equation for A is given by 

g^A = d'aeiZo, C, C) + ^^MZo, C, ()■ (4.4) 

It follows from the above equation that given any point .x" of the spacetime, the I = 2 and higher terms of 
Z are completely determined from asiu, C, C)- For example, up to linear order terms in as and/or cc' we 
have 



1 



(4.5) 

1 



= ^-f , (4.6) 

where the subscript R means the real part of the complex quantities. (If we keep bilinear terms of <tb and 
in 14.4! then the above terms also depend on a;'(u).) Inserting x*^ and x'^ in Eq.( l4.2l i gives an explicit 
relationship between and Re[ip^{u)Y. For the particular assumption given in ( 14.51 ). ( 14.61 ) we get 

It follows from the Bianchi identities that the time evolution of is quadratic in ctb. Assuming as, 
and vanish at = — oo and keeping up to second order terms in ctb or i?' we get the following 
expression 

D' + -^(j%P^ ^ MR\ (4.8) 
5v2c 

We have thus obtained an explicit relationship between the coordinate of the center of mass R^ and the Weyl 
scalars defined at Null Infinity. Even if we assume a more involved field equation for Z giving a functional 
dependence of .t'^ and x''^'^ on i?' as well as on as, Eq.( l4.2l i will give an algebraic expression relating i?* 
with the Weyl scalars at null infinity from which one can solve for the center of mass worldline. 

To obtain a relation between the velocity of the center of mass and the Bondi linear momentum we take 
a time derivative of ( I4.8l l. Using again the Bianchi identities yields the following expression 

+ r^^' " "''"''^^ " ^yitTn = "-^MV^ (4.9) 

(In the above expression we have omitted the terms A/, P', 5**^ , ^I''-''' since they are quadratic in ctb and (j)2.) 
Since all the vector quantities are aligned with the symmetry axis and the tensor variables are symmetric 
and trace free, we can get a simplified form of the equations as 

Z?^ + — ^cr|j^P^' Afi?^ (4.10) 

5v2c 

These equations provide explicit relations between D^, and R^, . 

The equation for the linear momentum can also be written as a sum of different parts as 



with 



= Pli+PB + PI,M, (4.12) 



Pli = (l-^a|f)^Ml^^ (4.13) 
" (4.14) 
Pmi = Y^'P14>IYr, (4.15) 



emphasizing the role of each contribution to the total linear momentum. The leading term in ( 14.151 ) is 
proportional to the charge times the time derivative of the electric dipole contribution. If one further assumes 
that the dipole contribution is due to a charged particle with worldline i?* one recovers a known result, the 
Abraham-Lorentz momentum0. However we are here concerned with astrophysical compact objects and 
this term usually vanishes. 

4.2 Equation of motion 

Taking a time derivative of ( 14.1 11 1 and inserting the Biachi identity 

= -^[7(<^'''^''+^''<^''') + ^(</'20§)1 (4.16) 

(4.17) 

gives the equation of motion for the center of mass, 

M{V''-^d^jlV'-) = -^[^{a^^^^i^-a=^¥''')^ + {o^-=^^=^)i,] (4.18) 

For completeness we also give the mass loss equation 

M = '^l^a^^w' + h^i'^w'^'] - (4.19) 

2V2G5 7 6c^^^^^ 

From the rh.s. of (14.18b we define the notion of gravitational and electromagnetic forces, i.e., 

Fg ^ ~^[lia^'''a^''~a^''^'"')n + i^^''^"')R] (4-20) 

Note that the gravitational force vanishes when the quadrupole or octupole moment vanish. In this case 
the electromagnetic radiation will produce the acceleration of the center of mass. This acceleration, how- 
ever, will be negligible for most situations. Since the mass loss equation has a separate contribution from the 
quadrupole and octupole moment, the net effect in this situation will be a reduction of the gravitational mass 
of the system while the center of mass remains at rest. On the other hand, most head-on collisions between 
compact objects will produce quadrupole and octupole radiation terms and there will be a net acceleration 
of the center of mass. 

It also follows from Eq. (14.18b that there are no runaway solutions. The functions afi, Fq and Fem 
decrease to zero as w oo. Asymptotically this equation gives a constant velocity if the total radiation is 
finite. Thus, the motion of the center of mass does not have runaway behavior 



5 Applications 

In this section we will first check that the formalism developed gives the correct answer for the cases 
where definitions have already been given. We will also compare our work with others to analyze simi- 
larities and differences in the definitions of linear and angular momentum. We begin with applying this 
formalism to the case of a stationary and axially symmetric spacetime. 

5.1 Stationary and axially symmetric spacetime 

First consider the Kerr metric, in this case we have V'4 = V-"" = 0- Moreover as the spacetime is stationary 
all derivatives respect ub vanish. In this way we would have to a = 0, so cr = (t(C, C)- Without loss of 
generality, we will work in a referential where cr = 0. From Eq. ( 12.37b we get that P' = vJ/'J = vj/^J'^' = 



for all i,j,k and from Eq. (12.311 ) we get that ip2 oc A/. Furthermore form Eq. (14.4b is easy to show that 

= (-V2R+t^a] (5.1) 



c2 V \/2c 

where a is the angular parameter The real part of ijj^^ is given from the Eq. (14.8b . so from Eq. ( IA.14b we 
have 

= aM (5.2) 

which corresponds to the angular momentum of the Kerr spacetime. If now we consider the Kerr-Newman 
case, the development is very similar to the Kerr case. Note that in Eq. ( IA.14b Acm'^ = 0, so we get the 
same equations 

= ^m(-V2R+i^^ (5.3) 
px ^ py ^pz (5 4) 

= aM (5.5) 

note that although these equations are identical, the evolution of the mass center is different for Kerr- 
Newman spacetime due to the presence of electromagnetic fields. 



5.2 Massive explosions or head on collisions. 

We consider here either a massive explosion of an isolated system, like type I supernova, or a massive 
head-on collision. We assume that initially the center of mass is at rest. Immediately after the explosion or 
collision the acceleration of the center of mass will be given by 

-^.[{mr + {mini 

Note that if either the quadrupole or octupole term vanishes there is no gravitational contribution to the 
acceleration. Any collision will have a quadrupole term but only collisions between uneven masses will 
also have an octupole contribution. Likewise, the electromagnetic force will be dominated by the radiation 
term since for most astrophysical objects = 0. 

Although total angular momentum is conserved, the coupling between gravitational and electromagnetic 
angular momentum gives a transfer mechanism by which the system can gain or loose intrinsic angular mo- 
mentum. Consider for simplicity that initially the system does not have angular momentum. After the 
explosion or collision the system will acquire an intrinsic gravitational angular momentum if electromag- 
netic radiation is emitted, i.e., from 

'^T = Jh-^W^iA^4>l)\i=i]=G- 

the electromagnetic angular momentum creates an intrinsic gravitational angular momentum in the opposite 
direction of the electromagnetic one. This effect could be important in charged isolated systems like the 
positron cloud discovered by the COMPTON detector in GRO, but will be negligible for most cases. 



5.3 Comparison with AKN equations 

In this subsection we will compare our equations with the Adamo-Kozameh-Newman (AKN) equations 
for linear and angular momentum. For simplicity we will consider a vacuum spacetime in the neighborhood 
of null infinity and assume the gravitational radiation only has quadrupole terms. Directly from [IJ we write 

P^ ^ ^McV^ ^^iV^-al,^ + V^-y^-) ~ ^^(QoY+Qaf) - l^(2i?V|j^ + ^laf ) (5.6) 



where the scalar ^| is related to the intrinsic angular momentum via 



Likewise, The total angular momentum in the AKN formalism is given by 

J| = ^McCI - - 2i?Vr) - ^MV^ar + ^^CI^F- (5.7) 

In this work the equivalent equations are given by 

P- = ^My-(l-^a^-) (5.8) 

and 

J^^S'' = const., (5.9) 

as one can see from applying the Komar approach and using a Bondi and a Newman-Unti cut as the bound- 
ary dTi. 

Although both formulations agree for stationary spacetimes, they differ when gravitational radiation is 
present. It follows from the Bianchi identities and the Komar integral for the axially symmetric Killing field 
that either in this or the AKN formulation we have 

P~ = const. 
J|, — const. 

Thus, in this formulation S^- ~ const, and V^- decreases to an asymptotic value after the gravitational 
radiation is emitted. In the AKN formulation both and S'^ are functions of time and obey a coupled 
system of ODEs. They also decay to an asymptotic value when the gravitational radiation is emitted. 

Also in this formulation the relationship between total and intrinsic angular momentum appears to be 
natural. Since the orbital part of the angular momentum vanishes when the center of mass vector i?* and the 
velocity V"^ are aligned along the z-axis, one expects that the intrinsic and total angular momentum should 
be equal. In the AKN formulation, the intrinsic angular momentum depends on the position and velocity 
of the center of mass and one does not expect this kind of relation. However, naturalness is not easy to 
define when gravitational radiation is present and it may well be that a more involved relationship of the 
form given in ( 15.7b is correct. A nice test for the two formulations will be available when gravitational wave 
astronomy is finally developed. 



6 Conclusions 

Using the available geometric structure of asymptotically flat spacetimes together with conservation laws 
that arise when those spacetimes are axially symmetric, we have defined the notion of linear and angular 
momentum for Einstein Maxwell spaces. 

Furthermore, using the light cone equation we have been able to identify worldlines inside the spacetime 
with Newman Unti cuts at null infinfity. The center of mass worldline i?° is then selected by imposing the 
condition that the mass dipole moment at null infinity vanishes when restricted to the center of mass NU 
cut. Using the available Bianchi identities at ,J^^ one obtains a relationship between the center of mass 
velocity and the Bondi momentum as well as the equation of motion of R'\ 

Several nice highlights of this approach are 

• A definition of angular momentum when electromagnetic fields are present. 

• A definition of center of mass worldline and velocity which are algebraically related to radiation 
fields at null infinity. 

• Definitions of gravitational and electromagnetic forces in terms of radiation fields. 

• Appropriate behaviour of the equations of motion (no runaway solutions). 



• A natural relationship between intrinsic and total angular momentum (they are the same in this case). 



The equations of motion could be used in astrophysical situations when the system has axial symmetry 
to predict the motion of the center of mass if the radiation is detected or to predict the amount of radiation 
if the velocity and acceleration of the center of mass is available. 

The formalism is ready to be generalized for spaces without symmetries and it will be considered in 
future work. In the generalization we expect some new features that are absent in axially symmetric spaces. 
Since at the moment there is no definition of angular momentum that has been universally accepted one can 
either work with a parameter dependent definition or use a suitable radiation data where all the definitions 
agree. It is left for the future to find a new definition of angular momentum for any kind of gravitational 
radiation. Acknowledgements: this research has been supported by grants from CONICET and the Agencia 
Nacional de Ciencia y Tecnologia. 



A Komar integral and angular momentum 

In vacuum spaces the 
angular momentum, 



In vacuum spaces the Komar integral of the Killing field ^^^^ yields a definition of the z-component of the 



J^ = -^lim (fv-ei^^dSab (A.l) 

St 

One can explicitly integrate this equation in the N-P formalism to obtain a formula at in terms of the 
spin coefficients. We first write the Killing vector field ^1*^^ as a combination of the null tetrad vectors as 

e'^) = 6^' + Crnm' + Urn'' + 6.^' (A.2) 

where 

= ^cos6' (A.3) 

r 

= ir COS 6 (A.4) 

C„ = (A.5) 

and the two-dimensional surface area can also be expressed as 

dSab ~ ^'^n^Jb]'!''^ sin OdOdip, (A. 6) 
thus the Komar integral can be written as 

J" = -— lim / / W''t^ JnJb~lanbysm0d(fde. (A.7) 

Using Eqs. ( 12.161 ), ( 12.201 ) and Eqs. ( I2.2l ),( l2.3l l and writing this equation up to order 0{r^^) we get 

cos 9 
r 

where w° = — 3(t° ll5l, lfT0ll . Thus, the Komar integral can be written as 



V^^C^) {nJb ~ lam) = ^/™(^? - frOcS") (A.8) 



1 

= - / Im{ij^ + a°Ba°)cos0d{cos0). (A.9) 
8 Jo 

where we have used the axial symmetry to integrate in the azimuth direction. Finally, this integral gives the 
following definition of angular momentum lfT3l 

J' oc Im{i'° + a°Ba°)\i=i (A. 10) 

We can follow a similar calculation with Eq. ( 13.81 ) 

= islToo/t^"^'-) + 2Ar^^,F^']dSab (A.11) 

St 



Using the fact that l5| 

F"'' = 2(/)oml"7i^l + + ml"™''!) + 2<?!)2;["m''l (A. 12) 

the second integral can be put in the form 

^^lini^ / 2A,Ci^^F''''dSab = ^^nn / 2A,^(^^),/)? smOdO (A.13) 



St 



so, we can define the total angular momentum as 



Jt = ^[Imi^i + + lim 1 2A,i? .4>\sinBdB] (A.14) 

12G r-+oo J 



Using the tetrad decomposition of the vector killing field ^^"^^ at ,J^+, and the fact that Ac and 0° are real, 
one can rewrite Jem as 

JiJA/ = lim / 2Ae^f^^(^?sin0d0 

= -lim / 2AX^ .^?d(cos0) 

= —lim / 2irAc(m'^ — m^)^5 cos6'd(cos6') 

= — 4 lim r / Im(Acm'^)(j)i cos 9d{cos 9) 
= -4 lim r/m(74cm''05)|/=i 



r^oo 

Furthermore, it can be shown that. 



lim rAcTn" = A" 

r— >-oo 



where is the free Maxwell potential data related to the electromagnetic radiation via 



00 = i°. 

Thus, the total angular momentum is finally expressed as 

B Tensorial spin-s harmonics products 

We present a table of tensorial harmonics products which complete the list of product |[l2 
Products of the form Yf^Y^J^l 



'-li ^3jkl 




2 p2(l) 
21 ^ijkl 


_ lzV2F.=^(V 
24 'J*:' 


_ ^^4(1) 

56 


(B.l) 


yOyl 


2\ vki 


8 p2(l) 
21 ^ijkl 


_ lzV2F.=^(V 
12 y*^' 


28 


(B.2) 


'^li-'Sjkl 


Y ijkl 


, 4 2(1) 


+ i*V2F.^(V ■ 


14 


(B.3) 




_ W 2(0) _ 
7 ijkl 


4^2(0) 


_ 1^^4(0) 
7 ^jkl 




(B.4) 


^li ^3]kl 


- A p2(0) _ 


2 ^2(0) 
21 ^^]kl 


+ lzV2F.t) 
24 


_ Ai;^4(0) 
56 


(B.5) 


^li-'Sjkl 


- Af2(o) _ 

21 ijkl 


Ar^^o) 


24^ y/^' 


_ Ai7^4(0) 

56 '^^^ 


(B.6) 



where 



ijkl 



ijkl 
ijkl 



^4ijkl 



Y-i^u„ + ^ilmY'^j]^^ 



with the superscript s = 0,1. 

Products of the form yslj^afc;™ 



yl -1^0 
^ 2ij ^ 3klm 

yO 

^2ij^?,klm 



^2ij ^3klm 



^ 2ij ^ 3klm 



^ 2ij ^ 3klm 



^ 2ij ^ 3klm 



^2ij ^3klm 



1 

12 



ijklm 



ijkhn 



1 



'l5 



.p3(l) 

ijklm 



2^ ijkhn 



' ^^^ijkhn ^'^^^ijklm + -^'^^^i\kh 



[^5(1) 
ijklm 



, /2^4(1) 2 

25 jjfeim ^ 24 'J*^'™ 21 

_ _ ^(-1(1) _ ^ H-i(i) + i^^/2F^*^^ - — iA/^G^^^^ 

-j^2 ijklm ijklm ijklm j ^ ijkhn ^ ijkhn 

+ , J_/^3(l) , J_7:t3(1) _ J_„-^/2f4(1) , J_p5(l) 

^ ijklm ' j^Q zj/c/m -j^Q ijklm 2g ijkhn ' -j^^ ijklm 

ijklm ijkhn ijklm 15S ijkhn ~^ 

_5_ /^„2(l) _J_ 1^3(1) _ J_(^3(l) „J_rr3(l) 

33g ^ ijkhn ^i*!:'™ 360 360 y'^'™ 

1680 u'fe'm 5040 'J*^'™ 

y ijklm ' '^i^^'^"' ijklm 

ijkhn 2Q ijklm 



14 



2g ijkhn 



30 



30 'J'^'^ 



ijklm 2^ ijklm -ij^/i 

+ ^iy3(i) +^iV2^^4(l) _^^5{l) 
gQ ijkhn 120 ijklm 210 '^3^^'^^ 

^ _ 1^(^1(0) _ 7^1(0) , ± p3(0) 

y ijklm ijkhn ijklm "ij/c/^ 

2 3(0) J_ nn5(0) 
15 '^jk^'f'^ 91 ijkhn 

2ij' 3khn y ijkhn ijklm ijkhn 

^„-^/2r2(0) , J_p3(0) 1 ^3(0) 

^ ijklm QQ-^iifelm "r X80 

LjV2F^^°^ _^F.5(o) 

168 u'*:'"! 252 'J*^'™ 

ijkhn ijklm ijkhn ^2 ijkh] 

J_p3(0) + J_(^3(0) 

'-"^ ijklm ' '^j^^' 



.^3(1) 

I ijklm 



1 

60 



;3(1) 

'ijklm 



15 



.(^3(0) 

ijklm 



!7n 



zlm 



180 



I ijklm 



7 



l7n 



— zx/2G'^°^ + — 
168 '^ ^■'"'^'^ 



zlm 



180 



.rr3(0) 

I ijklm 



, - . i7.5(0) 

A; / rri 252 '^3^^"^ 



^2ij ^3klm 



\r2 v — 2 



2^ ijkhn 2^ ijklm 



J_rrl(0) 

2^ ijklm 



1680 



ijklm ' 5040 



2^ ijklm 2^ ijklm 



2^ ijklm 



^ 336 ijklm 



1 



,3(0) 



1 



168 '•J''''™ 

^3(0) 1 



H. 



3(0) 



180 'J'^'™ 3go ^J*^'™ 350 ijklm 



1 



1680 



5(0) 



ijklm ^ 5040 



(B.15) 



(B.16) 



where (with ,s = 0, 1) 



l{s) 



G 

H. 



ijklm 
l{s) 
ijklm 

l(s) 



ijklm 
2{s) _ 



p 

ijklm 

(-<2(s) 

ijklm 



3(s 



G 
H 



ijklm 

3(s) 
ijklm 

3{s) 
ijklm 



F, 



ijklm 



F, 



5{s) 



[SikSjl + 5il5jk)yim + {^il^ jm. ^im ^jl 

(SjkSlm + SjlSmk + SjmSkl)Yii + {StkSlm + StlS„ik + StmSkl)Yij 
^vi^klY^m + SlmY^l + SmkY^i) 

{Slm£ikf + Skm^ilf + 5lk^imf)Y2jf + {Slm^jkf + Skm^jlf + 5lk^imf)Y2if 
Sili^jmfY^kf + <^jkfY2^j') + 5i,n{ejkfY2if + £]lfY2kf) + 
+ 6ik{^]lfY2mf + ^]mfY2if) + 5jl{timfY2kf + <^ikfY2mf) + 
+ 5jrn{^ikfY2if + eilfY2kf) + Sjki^ilfYi^f + eimfYiif) 

^ijYsklm 

^jmY-^ikl + ^jkY^iijn + ^jlY^ikm + ^imY'^jkl + ^ikY^jim + ^ilY^jkm, 
{^ikftjln + ejkf^iln + £ilfejkn + <^jlf<^ikn)Yr^f„^ + 

ilf^jmn + ^jlf^imn + ^imf^jln + ejmfeiln)Y^j:^f. + 
+ ( ^imf^jkn H" ^jmf^ikn H" ^ikf^jmn H" ^kmf^imnY^'^fnl 

^ikfyifjlrn + ^Hf^Ifjkm + ^im.fylfjkl + ^jkfylfUni + ^jlf^lfikm + ^j-mfYlfikl 



ijklm 



dijkhn 



Products of the form Y^^^j^Y^'^^^ 



yO yl 

3ijk 3lmn 



y-2 yi 
3ijk 3lmn 



rj^^^ijklmn 



y ^^^ijklmn 



-F, 



2(1) 



ijklm,n 



-21^- 
2 



2(1) 
ijklmn 



9. .^2(1) 
y ijklmn 



^'''^^ ijklmn 



15 ^4(1) 

154 '^3^^'^'^'^''- 



-^V2F'i'^ 

-j^g ijkLmn 

15 ^4(1) 
J^54 ijfc'm 



ijklmn 

25 4(1) 
154 ijklmn 



Liy^^5(l) Ai^6(l) 

V ~' ijklmn ' ]^^2 ijklmn 

—^^/2F^'^^^ - —iV2G^'^^^ +— f'^'^ + — G'^'^ 

140 ijklmn 2g ijklmn y2 ijklmn y2 ijklmn 



1 



-H. 



2(1) 



2(1) 

ijklmn 



252 y'^''"'" ' 56 

1 ;./2g3(l) _^p4(l) 1 p4(l) 

''TO ijklmn 515 ijklmn 5X5 ijklmn 

^^^V2F'!'l^+^F?^'^ 



—iV2F^i') - —iV2G''^y> 

215 tjkhnn 270 ijKimn 

5 



1848 



4(1) 
ijklmn 



6048 



ijklmn 15840 '^3^^'^' 



(B.17) 



(B.18) 



^3ijk-^3lmn 



yO yO 

3ijk 3lmn 



y-l T^l 



3ijk 3lmn 



+^^4(1) , 
2g ijkhnn ' ijklmn ' ^'j ijklmn ' 

^ ijklmn 528 '^3^^'^'^ 



14 

^2 ijklmn 

5 rr4(l) 
231 'i'jkhnn 



J_^2(l) 

ijklmn 



1 

14 



_i_H-2(l) 

42 y'^'" 



ijkhnn ^ ijklmn ' -j^2 ijkhnn ' -j^2 '^i^^'' 



2g ijkhnn 



1^^4(1) 

ijklmn 



-j^2 ijklmn 

J_(^4(l) 

ijklmn 



231 ijklmn 



1 



48 

"7 

8 

'21 
9 



ijklmn 

tr2(0) 

ijklmn 



■V2i^^'^' Lf6(i) 

^ 336 ijklmn 528 '^i^^^'''^ 

120^0 , lp2(0) ,1^2(0) 

y ^ijklmn ' ^ ijklmn ' ^ ijklmn 

y ijkhnn ijklmn 

__^4(0) 1^ rr4(0) ^ 6(0) 

ijklmn ijklmn 231 U^^^'^ 

rj^ijklmn'T' Y '-'ijklmn ^^''^ ijklmn ^ ^^'■^ '^'-'ijkhnn 

J_ p2(0) J_(-r2(0) _ J_ 0-2(0) A ^-2(0) 

-|^2 ijklmn -|^2 ijklmn 42 ijklmn 2g ijklmn 

'^^ijklmn ^ '^^^^ijklmn " '^^^^ijklmn 

1 p4(0) , 1 ^4(0) , 5 4(0) 
g]^g ijklmn ^ gj^g yfciimi ^ 1848 ijklmn 

^ ijklmn o/?r 



+ 



i 
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'3024 
2 



3696 



lO 



35 '^3^^'^^'^^ 
1 



<j ^ijklmn 



L J r^i 1 1 L 1 1 

'^^^^ijklm.n + Y4*"^^J's'i 



14 

'^^-^ijklmn ~ '^^^^^ijklmn ~ '^^'''^■^■^ ijklmn ~ 

J_^4(0) Lff4(0) 

880 ijkhnn 528 ijklmn 
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1 

880 Vklnm 



^ijklmn ~' 
ijklmn 
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7560 ijklmn 18480 y'^'™" 



ijklmn 
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l(s) 



ijklmn 
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ijklmn 



2{s) 



ijklmn 



G 
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ijklmn 
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Sil{SjmSkn + SjnSkm) + Sjl{SkmSin + Skn^im) + 5kl{5im5jn + 5in5j 

{SijSlm^knf + SijSln^kmf + SijSnm^klf + Sjk^lm^inf + SjkSln^imf H 
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+ {Sij5kn + SjkSin + SikSjn)Y2im 
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2(s) 
ijklmn 



ijklmn 



f: 



ijklmn 



3(s} 
ijklmn 



H, 



3{s) 
ijklmn 



ijklmn 

Mis) 

ijklmn 
ijklmn 



ijklmn 
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ijklmn 
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